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ABSTRACT 

We compute high temperature expansions of the 3-d Ising model using a recursive transfer- 
matrix algorithm and extend the expansion of the free energy to 24th order. Using ID-Pade 
and ratio methods, we extract the critical exponent of the specific heat to be a = 0.104(4). 
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1 INTRODUCTION 



High- and low-temperature expansions constitute major tools for the calculation of critical 
properties in statistical systems. The Ising and Potts model low temperature expansions 
were recently extended @, §, |, g using a technique based on the method of recursive 
counting ||. In a separate development, Vohwinkel Q implemented the shadow-lattice 
technique of Domb Q in a very clever way and added many new terms to the series. 
However, the extraction of critical parameters from low temperature series is hampered 
by the presence of unphysical singularities. This is especially true of the 3-d Ising model. 
For this reason, low temperature analytic methods are very often inferior to Monte-Carlo 
methods for computing critical exponents. 

High-temperature (HT) expansions on the other hand, generally have better analytic 
behavior and yield more accurate exponents. Very recently, two variants of the recursive 
counting technique for HT expansions have been pursued. While Enting and Guttmann 0] 
keep track of spin configurations on a set of rectangular finite lattices, ref-Q counts HT- 
graphs on finite, helical lattices. Such computer based series expansions have very large 
memory requirements. This makes them ideal candidates for large parallel computers if 
communication issues can be handled efficiently. In this paper we will present the results 
of a HT expansion of the 3-d Ising model to 24th order, obtained on a 32 node 1 GByte 
Connection Machine CM-5. The implementation is based on a bookkeeping algorithm of 
binary coded spin configurations in helical geometry. 



2 COMPUTATION OF THE SERIES 

We start with a discussion of the HT algorithm to compute the partition functions on finite 
3-d Ising lattices. Starting from the action 

E{s} = -J2 s i s J > (!) 
<i,j> 



the partition function is 



Z = exp (-(3E) = ]T J] ex P (P*i°j) ( 2 ) 

{s} {s} <i,j> 

and is expanded in a HT series [|| 

Z = (cosh 3 (3) v J2 II ( X + **J*) = ( 2 cosh3 , (3) 

{s} <i,j> k 

with the HT expansion parameter t = tanh (3, V is the volume of the system. The free 
energy per spin is defined as 

/ = -^iogz = - 2cosh3/3 -ly / ^. ( 4) 

/3V p 

For simplicity, consider a finite simple cubic lattice which, in the recursion algorithm, is 
built up by adding one site after the other, layer by layer. This procedure defines the 
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recursion step, which requires knowledge only of those spin states that are contained in 
the exposed two-dimensional surface layer. To minimize finite size effects, it is best to 
use helical boundary conditions [|. One can visualize helical boundary conditions by 
imagining all spins in the layer laid out along a straight line. In this picture, the nearest 
neighbours to a given site in the sequence in the ith direction can be chosen to be hi sites 
away, with i = x,y,z. It is convenient to assume h x < h y < h z . It is easy to see that as 
spins are added, one needs only to keep track of the states of spins on the topmost h z sites. 
Let these spins be denoted st, . . . , Sh z . Then the partition function can be rewritten as 



Z = (2cosh 3 /3) y ]T P(k;s 1 ,...,s hz )t h 

k s 1 ,...,s h 



(5) 



The recursion step, which consists of adding another spin so to the system, changes the 
partition function into 



Z = 2 v '(cosh 3 /?) v '+ 1 ££ £ p(k; Sl ,...,s h Jt k 

s Q k s 1 ,...,s hz 
X (1 + SQS hx t)(l + S S hy t)(l + S S hz t) . 



(6) 



Sfi x , Sh y and Sh z are the backward nearest neighbours of the site sq. The site so will displace 
its backward z neighbour site Sh z after the counting of the added spin is completed. Since 
Sh z will not be referred to in the subsequent steps of the algorithm, the summation over 
Sh z can be carried out: 



Z = 2^(cosh 3 /3)^£E £ 

s k s 1 ,...,s hz _ 1 

x [ +p(k; si, . . . , s^_i, s )t k (l + s s hx t)(l + s s hy t)(l + t) 
+p(k; st,..., s hz - 1 ,s )t k (l + s s hx t)(l + s s hy t)(l - t)] . 



(7) 



The contribution in the second (third) line of this equation contains the part with Sh z being 
parallel (antiparallel, denoted by s~o) to sq. Comparing this expression with the HT series 
@ for the new system yields the recursion relation induced for the coefficients p: 

2p'(k;s ,st, ■ ■ ■ ,s hz -t) 
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It is crucial to remove finite-size errors by combining the results of different lattice structures 
as described in refs. |2|, ||. We use the set of lattices listed in table |TJ and obtain the 
free energy coefficients up to 24th order as given in table |2[ In order to eliminate the 
contribution from (unphysical) loops with an odd number of links in any direction, we use 
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Table 1: Structures and weights w of the lattices used 
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Table 2: Free energy up to 24th order 
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the cancellation technique of ref. || . This amounts to inserting additional signature factors 
into eq. || for each of the three link-factors 

(1 + s s h j) -> crj(l + s s hx t) , i = x,y,z (9) 

with {cr x ,(Jy, cr z } = {±, ±, ±}. By performing 8 separate runs corresponding to all possible 
values of a and adding the results, one achieves a complete elimination of the unwanted 
loops. Possible contributions of higher-order finite-size-loops are at least of order 25 for this 
set of lattices. Since we use open boundary conditions, the coefficients p are invariant under 
the global transformation Sj — > — Sj. This Z(2) symmetry enables us to reduce memory 
requirements by a factor of two. Unlike refs. [||, ||, ||] we use multiple-word arithmetic 
to account for the size of the coefficients. This implementation needs about 100% more 
memory but leads to a doubling in performance. Since the number of words can be adjusted 
separately for every order, the computational effort can be reduced accordingly. On the 32 
node CM-5 the total time for all computations was about 50 hours. 

Compared to the finite-lattice approach of Enting and Guttmann M, our method ap- 
pears to require more CPU-time since we need to cancel unphysical loops. It should be 
noted, however, that helical lattices are very naturally implemented in data parallel soft- 
ware environments and thus lead to better performance. In the usual finite lattice method 

the HT expansion can only be extended in fairly coarse steps, using lattices with (4 x 5) 
cross-section for 22nd order and (5 x 5) cross-section for 26th order, respectively. For this 
reason, a 24th order computation would not have been feasible using that method with our 
computer resources. 



3 CRITICAL EXPONENT 

The specific heat is defined as 



d 2 



c\h=o = P 2 g^ l °S z = J2 c kt 



2k 



and is expected to behave near Tq as 

c \ h=0 = A(T)\T-T c \ 



1 + B(T)\T -T c \ e + ... 



(10) 



(11) 



with A and B being analytic near T c ||, 0. We analyse the series using unbiased and 
biased inhomogeneous differential Pade-approximants (IDPs) fill] as well as ratio-tests. 



3.1 Pade- Analysis 

In figure |l| we plot a against t^. for each IDP-Approximant [J/L;M]. Fitting the linear 
dependence of a on Q, we find 

a = 0.102 ±0.008 (12) 

at the value tc = 0.218092 as obtained in Monte-Carlo simulations [jT^ . A direct, biased- 
IDP analysis was also performed. We obtained a = 0.109 ± 0.016. 
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Figure 1: Critical exponent a as a function of t^. 
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IDPs can also be used to predict the most significant digits of the next term in the 
specific heat series Q. The estimate of the 24th order term as obtained in ref. [||| agrees 
perfectly with our exact result. Using the same method we can estimate the 26th order 
term in the expansion to be 

/ 26 = 443762(4) x 10 7 , (13) 
where the errors quoted are two standard deviations. 

3.2 Ratio-Test 

The main problem in the determination of critical exponents in the low-temperature case 
is the presence of unphysical singularities nearer to the origin than the physical one. Since 
the expansion coefficients c n are dominated by these unphysical singularities, ratio-methods 
cannot be applied. 

In the HT-expansion, the physical singularity dominates the asymptotic behaviour, so 
that the ratio r n = c n /c n _i of successive coefficients of the series is expected to behave as 



n ~ 4.2 



If a - 1 c d n ( 1 



n 



n 



1+26 



n l+36» 



(14) 



Assuming that the correction-to-scaling exponent 9 is close to 0.5 12, 13 1, the following 
sequence s n is expected to converge towards a like 



3 n ; — (tc r n ~ l) 



l)n+l = a + 



c d / 1 

n 1 / 2 n \n 3 / 2 



(15) 



A plot of this sequence against n is shown in figure |2[ Obviously the first four values are 
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Figure 2: Plot of the sequence s n against n. The error of S13, obtained from the ID-Pade 
extrapolation, is too small to be visible. 

dominated by higher order corrections. To obtain estimates for a we therefore use only 
the values {sq, . . . , S13}. A 3-parameter least-square-fit using the ansatz of eq. fli~5| ) yields 
the values shown as diamonds in figure [3|. The value of a = 0.113 obtained by the fit to 
the points {sq, . . . , sn} is in perfect agreement with the result of ref. Their estimate 
of a = 0.110 using the extrapolated term s\2 appears to be slightly above our value of 
a = 0.108 using the exact term. Including our value for S13 of the ID-Pade extrapolation 
eq. ([i~3| ) we obtain a = 0.105(2). The error represents the uncertainty of the extrapolation. 
However, from fig. ^ it it quite suggestive that the a- values might converge to a value 
below 0.105. 

To get an estimate of the uncertainties of our results, we investigate the stability of the 
fits. For this purpose, we repeat the analysis after eliminating the point sq from the data. 
As a result we obtain sizeable changes for a. The new data are shown as crosses in fig. |3| 

In figure |] we present the results for the first correction-to-scaling coefficient c from 
our 3-parameter fits. In contrast to ref. Jl, our values suggest that c changes sign with 
increasing n max . Because of the sensitivity of the fits to the number of terms we keep, it 
is difficult to determine the value of c very precisely. Our best estimate is c = 0.01(4). 
Since c vanishes within error, it seems reasonable to also try a 2-parameter-ansatz with 
c = to fit the data. The results of these fits are shown in figure |. We now find that the 
fits are much more stable and the a estimates show much more of a convergence to their 
asymptotic values. The best value (from the largest n max ) is a = 0.1045(3). This value 
supports the impression of the 3-parameter fits, which suggested that a was slightly below 
0.105. Taking into account the fact that neglecting c causes a systematic error, our final 
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Figure 3: Estimates of alpha using a 3-parameter-fit. Each point represents the results of a 
fit to the set of values {s nmin , ■ ■ ■ , Sn max }- The error bars of the rightmost values represent 
the uncertainty of the extrapolated 13th term. 
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Figure 4: Estimates of c using a 3-parameter-fit. Each point represents the results of a fit 
to the set of values {s nmin , • • • , s nmax }. The error bars of the rightmost values represent the 
uncertainty of the extrapolated 13th term. 
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Figure 5: Estimates of a using a 2-parameter-fit with c = 0. Each point represents the 
results of a fit to the set of values {s nmin , . . . , s nmax }. The error bars of the rightmost values 
represent the uncertainty of the extrapolated 13th term. 



estimate for the critical exponent is, 

a = 0.104(4) . (16) 

4 DISCUSSION AND OUTLOOK 

The crucial element in the estimate of the error in a ( eq. |l6| ) is our neglect of the 
correction-to-scaling coefficient c. The resulting systematic error is rather large. From fig. 
H one might speculate that the estimates for c begin to exhibit asymptotic behaviour at the 
26th order. Therefore an exact calculation of the 26-th term of the expansion might reduce 
the uncertainty of c significantly. If the magnitude of c turns out to be really negligible, one 
could adopt the errors of the linear fits, and a would be obtained accurate to the fourth 
significant digit. 
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